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The electromagnetic form factors of the spin-3/2 D baryons, namely D, fl/, it*,, and Qccc, are 
calculated in full QCD on 32® x 64 PACS-CS lattices with a pion mass of 156(9) MeV. The electric 
charge radii and magnetic moments from the EO and Ml multipole form factors are extracted. 
Results for the electric quadrupole form factors, E2, are also given. Quark sector contributions 
are computed individually for each observable and then combined to obtain the baryon properties. 
We find that the charm quark contributions are systematically smaller than the strange-quark 
contributions in the case of the charge radii and magnetic moments. E2 moments of the 12/^ and 
flccc provide a statistically significant data to conclude that their electric charge distributions are 
deformed to an oblate shape. Properties of the spin-1/2 flc and Dec baryons are also computed and 
a thorough comparison is given. This complete study gives valuable hints about the heavy-quark 
dynamics in charmed hadrons. 

PACS numbers: 14.20.Lq, 14.20.Jn, 12.38.Gc, 13.40.Gp 


I. INTRODUCTION 

Since the discovery of the proton’s internal structure, 
form factors have been the tools to investigate the in¬ 
ner structure of hadrons. One of the intriguing proper¬ 
ties of hadrons is their electromagnetic structure, such 
as their charge radii and magnetic moments. There have 
been enormous efforts to determine the electromagnetic 
form factors both experimentally and theoretically. On 
the theoretical side, lattice QCD is a widely used first- 
principles calculation framework to study these form fac¬ 
tors. Main challenges for the lattice QCD form factor 
calculations have been the pseudoscalar/vector-meson 
states and the nucleon (see [I] for a review) with only 
recent works probing the octet [2-4] and decuplet sec¬ 
tors [5-7]. 

Improvement in our understanding of the light hadron 
electromagnetic structure makes it timely to study the 
heavy-flavor hadrons further. Comparison of the two 
sectors would reveal differences in the quark-gluon dy¬ 
namics of heavy flavors. We have recently examined the 
internal structure and the quark dynamics of the hadrons 
that contain at least one charm quark using lattice QCD. 
Initially in the meson sector, we found that the electric 
charge radii and the magnetic moments of the D and 
D* mesons to be smaller [8] as compared to those of the 
TT and p mesons. An investigation of the quark contri¬ 
butions revealed that the decrease is mainly due to the 
smallness of the charm-quark contributions to the ob¬ 
servables. In this same vein, we extended our calcula¬ 
tions to the spin-1/2 light-charmed and strange-charmed 
baryons. We observed a similar behavior in the inter¬ 
nal structure caused by the heavy charm quark [9, 10]. 
Finally, we made a lattice study of the experimentally 
observed Q* t ( 5 ^) 7 radiative transition and 


estimated its lifetime [11]. 

In this work, our aim is to broaden our perspective 
by including the elastic electromagnetic form factors of 
the </ = f strange-charmed baryons to examine the 
spin-dependence of the quark dynamics. We compute the 
electromagnetic form factors of the Q, Q*, Q/c and Occe 
baryons as well as the Oc; ^cc baryons with J = . 

We extract the electric charge radii and the magnetic 
moments, and give a thorough comparison of both for 
the spin-3/2 and spin-1/2 sectors, which helps to improve 
our understanding of the nonperturbative structure of 
the strange-charmed baryons. 

Our work is organised as follows: In Sec. II we outline 
the continuum and the lattice formulation to calculate 
the electromagnetic form factors of the spin-3/2 baryons 
and discuss the lattice setup and the methods we utilise. 
Sec. Ill is devoted to the results. We present our results 
for the baryon masses, electric charge radii and magnetic 
moments of the spin-3/2 baryons as well as the results 
for their spin-1/2 counterparts. The electric quadrupole 
moments of the spin-3/2 baryons are also given in this 
section. Our main findings are itemised at the end of 
each subsection and a summary is given in Sec. IV. 


II. THEORETICAL FORMULATION AND 
LATTICE SETUP 

Electromagnetic form factors of baryons can be calcu¬ 
lated through their matrix elements of the electromag¬ 
netic vector current ^ eqq{x)^^q{x), where q runs 
over the quark content of the baryon in consideration. 
We refer the reader to Ref. [10] for the analytical expres¬ 
sions of spin-1/2 baryons. Here we only give the details 
for spin-3/2 baryons. 
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A. Theoretical formulation 

The electromagnetic transition matrix element for the 
spin-3/2 baryons can be written as 

s)} 

I m2 (1) 

= Y ,s')0'^>^^Ur{p,s), 

where p{s) and p'{s') denote the four momentum (spin) 
of the initial and hnal states, respectively. Mg is the 
mass of the baryon, E {E') is the energy of the incoming 
(outgoing) baryon state and Mq(p, s) is the baryon spinor 
in the Rarita-Schwinger formalism. The tensor in Eq. (1) 
can be written in a Lorentz-covariant form as [12] 




(2Mb)2 




where P = p+p' and q = p' —p. The multipole form fac¬ 
tors are defined in terms of the covariant vertex functions 
oi, 02 , Cl and C 2 as, 

Geo{Q^) = (1 + o ''") Wl + (1 +t)0'2} 

1 ^ (3) 

- 3'r(l + 'r){ci-b(l-fT)c2}, 

GE2{q^) = {ai + (1 + 'r)o2} 

1 (4) 

- 2^^ + r) {ci + {1 + t)c2} , 

GMi(g^) = (1 + ^T)ai - ^t( 1-bT)ci, (5) 


where at are the Pauli spin matrices, a, /3 denote the 
Dirac indices and cr, r are the Lorentz indices of the 
spin-3/2 interpolating fields. The baryon interpolating 
fields are chosen, similarly to those of Delta baryon as 

= ^P^’^{2[ql\x)C-/^(i{x)]q^{x) 

v3 (10) 

+ [ql\x)Cjf,qi{x)]ql^{x)}, 

where i, j, k denote the color indices and C = 
7472 - <?i, <? 2 , <?3 are the quark flavors and chosen as 

( 91 , 92 , < 73 ) ={(s,s,s), (s,s,c), (s,c,c), (c,c,c)} for O, D/, 
and flccc baryons, respectively. It has been shown in 
Refs. [5, 13] that the interpolating field in Eq. (10) has 
minimal overlap with spin-1/2 states and therefore spin- 
3/2 projection is not necessary. 

Inserting a complete set of eigenstates 

|(p, s))((p, s)| into (7) and (8) and taking the 
large Euclidean time limit, t 2 — ti and ti ^ a, 
correlation functions reduce to 

(G»f(t;p;r4)) = 

ZBip)ZE{p)^e-^^TT[T,A,,], 

][I‘2 

(G®f®(t2,ti;p',p;r)) = Zb{p')Zb{p)^ 


{'^‘2.—tl) ^ — Et\ 


Tr[rA,,,(p')G"''" A,,,(p)], 


where the trace acts in the Dirac space, the Zb{p) is 
the overlap factor of the interpolating field to the corre¬ 
sponding baryon state and A^t is the Rarita-Schwinger 
spin sum for the spin-3/2 field in Euclidean space, defined 


= ai--{1 +t)ci, (6) 

with T = —(7^/(2 Mb)^. These multipole form factors are 
referred to as electric-charge (£10), electric-quadrupole 
{E2), magnetic-dipole (Ml) and magnetic-octupole (M3) 
multipole form factors. 

The two- and three-point correlation functions for spin- 
3/2 baryons are defined as, 

(G«»(t;p;r4))=^e-*PMr' 

X ^ (7) 

X (vaclT[r7“(a;)7?“'(0)]lvac), 


(G»f«(t 2 ,ti;p',p;r)) = -z ^ 

X2,Xi 

X r““' (vaclr[r7“ ( 0:2 )jf, {xi)fi“' (0)] [vac), 
with the spin projection matrices 


0\ 

* 2 1 0 or 


'’• = 5(0 o)- 


'^Ua{p,s)Ur{p,s) = 


—17 • p + Mb 


X ygar + 3^2 » 

= A,j7.(p). 


2Mb 

^PaPr _ . Pair - PtIA ( 13 ) 

3M| * 3Mb 


To extract the multipole form factors we consider 
the following ratio of the correlation functions given in 
Eqs. (7) and (8), 


i?^^r(i2,ti;p',p;r) = 

- p. r))(Gg/'‘^(f 2 , G; p, -p'; r)) 1 

(GfB(t2;p';r4))(Gfe(t2;-p;r4)) 
ti»a , (E^ + MbV'^ (Ep,+MbV'^ 


t2-ti»o V 2E, 


xn,^(p',p;r). 


Note that there is no sum over the repeated indices. 
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The multipole form factors can be extracted by using 
the following combinations of 11^p; F) [5]: 


G£o(<z') = I (nA(q.,0;r4) + ^A(q^,0;^4) 
+ 113^3 (qi, 0 ; r 4 )), 


(15) 


GE2{q^) = 2 


2 , MiE + M),^ 


(nA(q.,0;r4) 


+112 2(qi; Oj ^4) 2113 3(qi, 0; r4) ), 


(16) 


GMiiq^) = (ni\(qi,0;r2) 

+n2^2(qi> 0; + 2 ) + n3^3(qi, O; r2)), 


(17) 


Gmsiq^) — “4 


M{E + Mf 


|qil 


(n3\(qi,0;r2) 


(18) 


+n2^2(qi) 0; + 2 ) — 2n3^3(qi> O; + 2 )), 


where i = 1, 2, 3 and qi are the momentum vectors in 
three spatial directions. In case of the E2 form factor, it 
is possible to exploit the symmetry, 

n242(q„0;r4) = n343(q„0;r4), (19) 


and define an average 
nLg(q*,0;r4) = 

\ [n2'‘2(q*,0;r4)+ n3‘‘3(qi,0;r4)], 


( 20 ) 


in order to decrease the statistical noise in Ge 2 - With 
the above definitions, Ge 2 form factor can be rewritten 
as 


Ge2{q ) = 2 


2 , ,M(i? + M). 4 


-n^ 

^-^avg 


|q^r 
(q,,0;r4)) . 


(^A(q^,0;^4) 


( 21 ) 


We consider an average over momentum directions for 
both EO and E2 form factors. In case of the Ml form 
factor, we make a redefinition to utilise all possible index 
combinations in order to improve the signal. Sum of all 
correlation-function ratios for Ml is written as 


Gmi(<7^) 


3{E + M)1 
5 |qp 6 


3 

E 


n4,(qi,0;rfe) 


i^j,k=l 


( 22 ) 


+niVqfcj0;O + ni*i(qj;0;rfe) 


Compared to the dominant form factors EO and Ml 
we have observed that the data for the E2 and M3 form 
factor is much noisier. It turns out that with the limited 
number of gauge configurations we have at the smallest 


quark mass, the data for M3 moments are too noisy to 
allow a statistically significant value. Thus, we omit the 
M3 form factor in this work and extract only the EO, 
Ml and E2 form factors for the lowest allowed lattice 
momentum transfer. 

It is possible that the higher order form factors in the 
expansion interfere with the leading and sub-leading form 
factors that we consider. Although a dedicated study 
would be needed to have a strong conclusion we note 
that the agreement between the results obtained from 
different lattice formulations of Ref. [5] and Ref. [7] sug¬ 
gests that the interference effects are minimal. 


B. Lattice setup 

We have run our simulations on gauge configurations 
generated by PACS-CS collaboration [14] with the non- 
perturbatively 0(a)-improved Wilson quark action and 
the Iwasaki gauge action. The details of the gauge config¬ 
urations are given in Table I. The simulations are carried 
out with near physical u,d sea quarks of hopping param¬ 
eter = 0.13781. This corresponds to a pion mass 
of approximately 156 MeV [14]. The hopping parameter 
for the sea s quark is fixed to = 0.13640 and the 
hopping parameter for the valence s-quark is taken to be 
the same. 


TABLE I. The details of the gauge configurations used in this 
work [14]. We list the number of flavors (Nf), the lattice spac¬ 
ing (a), the lattice size (L), inverse gauge coupling (/3), clover 
coefficient (csvv), number of gauge configurations employed 
(Ngc) and the corresponding pion mass (m,r). 


Nl X Nt Nf a [fm] L [fm] 

P 

Csw Ngc [MeV] 

32® X 64 2+1 0.0907(13) 2.90 

1.90 

1.715 194 156(7)(2) 


For the charm quarks, we employ the Clover action and 
use the hopping parameter value, Kc = 0.1246, which 
we have determined in our previous work [10]. In or¬ 
der to tune the hopping parameter we apply the Fer- 
milab method [15] in the form employed by the Fermi- 
lab Lattice and MILC Collaborations [16, 17]. A similar 
procedure has been recently used to study charmonium, 
heavy-light meson resonances and their scattering with 
pion and kaon [18-20]. In the Fermilab method’s sim¬ 
plest application one sets the Clover coefficients ce and 
cb to the tadpole-improved value I/mq, where uq is the 
average link. We follow the approach used in Ref. [18] to 
estimate the uq as the fourth root of the average plaque- 
tte and determine the charm-quark hopping parameter 
Kc nonperturbatively by tuning the spin-averaged static 
masses of charmonium and heavy-light mesons to their 
experimental results. 
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We make our simulations with the lowest allowed lat¬ 
tice momentum transfer q = 2t:/[N go), where Ng is the 
spatial dimension of the lattice and a is the lattice spac¬ 
ing. This corresponds to three-momentum squared value 
of = 0.183 GeV^. In order to increase statistics, 
we insert all possible momentum components, namely 
{\qAMv\Mz\) = (-1,0,0), (0,-l,0), (0,0,-l), (1,0,0), 
(0,1,0), (0,0,1). We also consider vector-current and 
spin projections along all spatial directions and take into 
account all Lorentz components of the Rarita-Schwinger 
field. We employ a wall-source/sink method [8], which 
enables us to simultaneously extract all the components 
of the correlators given in Eqs. (15-17), eliminating the 
need for extra inversions. However, since wall smearing 
is not a gauge-invariant smearing method, gauge hxing 
becomes a necessity. We choose to work with Coulomb 
gauged configurations, which lead to a better ground- 
state saturation. 

In order to improve the ground-state coupling, non¬ 
wall smeared source and sink are smeared in a gauge- 
invariant manner using a Gaussian form. In the case of s 
quark, we choose the smearing parameters so as to give 
a root-mean-square radius of (n) ~ 0.5 fm. We have 
measured the size of the charm quark charge radius to 
be small compared to the light and strange quarks, both 
in mesons [8] and baryons [9, 10]. Therefore, we adjust 
the smearing parameters to obtain (rc) = (r;)/3. 

The source-sink time separation is fixed to 1.09 fm 
(t 2 = 12a), which has been shown to be sufficient to 
avoid excited state contaminations for electromagnetic 
form factors [10]. Using translational symmetry, we have 
employed multiple source-sink pairs by shifting them 12 
lattice units in the temporal direction. All statistical er¬ 
rors are estimated by a single-elimination jackknife anal¬ 
ysis. We calculate the connected diagrams only and con¬ 
sider the point-split lattice current 

(23) 

which is conserved by Wilson fermions. 

III. RESULTS AND DISCUSSION 
A. Baryon masses 

The masses of the U, flc , Ucc and flccc baryons are 
extracted from the shell-source/point-sink lattice two- 
point function given in Eq.(7) by a simultaneous fit to 
all spatial Lorentz components. We show the effective 
mass plots for the spin-3/2 baryons in Fig 1. Results are 
given in Table II along with a comparison to the masses 
reported by PDG and other lattice collaborations. 

We identify the lower end, tmm, of the fit regions 
[tmin,tmax] by Searching for a time slice where a plateau 
forms in the effective mass plots. The upper end, tmaa,, 
extends up to time slices where the signal starts to dete¬ 
riorate. We note that shifting tmin to larger time slices 


gives lower mass values that are closer to those of other 
lattice groups or the experimental ones. However, the ef¬ 
fective mass plots are known to exhibit a fluctuation after 
a plateau is formed [21] and shifting the fit region would 
introduce a bias to the determination of the masses. 

As compared to the experimentally available results 
there is a discrepancy of around 100 MeV in the case of 
U and flc masses. Note that the differences may arise 
from our choice of the strange and charm quark hopping 
parameters. In order to avoid the partial-quenching ef¬ 
fects we have chosen Kg to be the same as that of the 
sea quark. On the other hand, our 0 mass is in good 
agreement with the mass reported by the PACS-CS Col¬ 
laboration [14]. A retuning of Kg so as to obtain the 
physical K mass would be desirable for precision calcula¬ 
tions. However we expect such a retuning to have a min¬ 
imal effect on the conclusions of this work. In Ref. [9], 
by re-tuning Kc mildly, we have indeed confirmed that a 
variation of baryon masses around ±100 MeV has mini¬ 
mal effect on the form factors. We have found that such 
a tuning changes the charge radii by less than 2%. 

The mass of the VLccc as obtained in our simulations 
agrees with those from other lattice simulations [13, 22- 
24] with different actions. This can be taken as a good 
indicator for the aptness of the charm-quark action we 
employ. 



FIG. 1. Effective mass plots for the spin-3/2 baryons. Dashed 
horizontal lines indicate the fit regions. 


B. Form Factor correlation functions 

In our simulations we evaluate each quark sector sep¬ 
arately and normalise to unit charge contributions. The 
baryon properties are then estimated by combining the 
quark contributions with their weights from respective 
quark numbers and electric charges as, 

{0) = Ngeg{0g)+N,e,{0^), 


(24) 
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TABLE II. Q, n, (| + , f+), n,, (| + , 1+) and flccc masses (at a pion mass of m,r = 156 MeV) together with the experimental 
values [25] and those obtained by PACS-CS [22] (at the physical point, except 11 which is at m,r = 156 MeV [14]). We have 
also included results by ETMC ]13], Briceno et al. [23] and Brown et al. ]24]. All values are given in units of GeV. 




This work 

PAGS-GS [22] 

ETMG [13] 

Briceno et al. [23] 

Brown et al. [24] 

Exp. ]25] 



[GeV] 

[GeV] 

[GeV] 

[GeV] 

[GeV] 


ric 

1 + 

2 

2.783(13) 

2.673(17) 

2.629(22) 

2.681(48) 

2.679(57) 

2.695(2) 


1 + 

2 

3.747(10) 

3.704(21) 

3.654(18) 

3.679(62) 

3.738(40) 

— 

Q. 

3 + 

2 

1.790(17) 

1.772(7) [14] 

1.672(18) 

— 

— 

1.673(29) 


3 + 

2 

2.837(18) 

2.738(17) 

2.709(26) 

2.764(49) 

2.755(61) 

2.766(2) 


3 + 

2 

3.819(10) 

3.779(23) 

3.724(21) 

3.765(65) 

3.822(42) 

— 

f^cce 

3 + 

2 

4.769(6) 

4.789(27) 

4.733(18) 

4.761(79) 

4.796(26) 

— 


where (O) is the observable, Nq is the number quarks 
inside the baryon having flavor q and Cq is the electric 
charge of the quark. 

We have extracted the spin-3/2 baryon multipole form 
factor values by searching for plateau regions of the ra¬ 
tio given in Eq.(14). The correlation-function ratios for 
the EO, Ml and E2 form factors are depicted in Figs. 2- 
4. Fit values for the form factors at the lowest allowed 
three-momentum transfer (q^=0.183 GeV^) are given in 
Table III. Note that EQ form factor reduces to the elec¬ 
tric charge of the baryon as usual and the other form 
factors cannot be directly obtained at zero momentum 
transfer due to their definitions in Fqs. (16)-(18). 


C. Charge radii 


Electric charge radius of the baryons are obtained 
by calculating the slope of the EO form factor at zero- 
momentum transfer: 

{r%) = -Q^GEom\Q^=o- (25) 

In the case of the proton, the low-Q^ experimental data 
is well-described by the dipole form Ansatz 

= (1 (20 

where A is the dipole mass. We assume that such Ansatz 
also holds for the baryons we study here. Since we per¬ 
form our simulations with a single value of the finite mo¬ 
mentum transfer, a dipole fit of the form factor to a mo¬ 
mentum region is not possible. We can, however, extract 
the charge radii using the expression 


( 4 ) 12 ( I Geo{0) 

Geo{ 0) QLn \\ GEoiQLn) 


( 27 ) 


which can be readily derived by inserting Eq.(26) into 
Eq.(25). Our numerical values for the electric charge 
radii are given in Table IV. Note that the quark sector 
contributions are for individual quarks of unit electric 
charge. 

We observe that the s-quark contribution to the elec¬ 
tric charge radii is almost independent of the quark-flavor 
composition of the baryon. The charge radii of both spin- 
1/2 and spin-3/2 baryons agree within one standard de¬ 
viation, which can be seen more clearly in Fig. 5. In the 
case of c-quark contributions illustrated in Fig. 6, the 
charge radii in all baryons are similar. Although the con¬ 
tribution of individual c-quark seems to increase slightly 
as the number of c-quark increases in the composition of 
the baryon, this change remains negligibly small. 

In Fig. 7 we show the ratios of individual quark-flavor 
contributions in the spin-1/2 to that in the spin-3/2 sec¬ 
tor. We observe that for the singly-charmed baryon, 
s- and c-quark charge distributions are insensitive to the 
spin-flip of the c-quark whereas in the case of the doubly- 
charmed rice baryon the contributions of s- and c-quark 
to the charge radii increase. 

We combine the individual quark contributions accord¬ 
ing to Eq. (24) and list the numerical results in Table IV. 
We find the electric charge radius of the fl baryon to be 
(r|;)f 2 - = —0.326(21) fm^ in quite good agreement with 
the previous lattice determinations [5, 7]. A comparison 
of baryon charge radii is made in Fig. 8. In magnitude, 
n baryon has the largest electric charge radius among 
all baryons we study. Spin-1/2 (spin-3/2) Qc (f2*) and 
riccc seem to have similar charge radii while the Qcc (S^cc) 
baryon has almost a vanishing charge radius. 

Based on the similarity of the quark contributions 
to the charge radii one can naively assume that the 
quark sector contributions to the charge radii to be sim¬ 
ilar for all spin-3/2 baryons that we consider so that, 

(^l)n = = 

= ^c- Using Eq. (24) we can derive a rela- 
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TABLE III. The values of E0{Q^), M1{Q^) and E2{Q^) form factors at q^=0.183 GeV^ for fl, n*, and Qccc- Results are 
given in lattice units for single quark and normalised to unit charge. 



A0"(Q") 

A0"(Q") 

MViQ'^) 

MR(Q") 

E2‘’iQ'^) 

E2%Q'^) 

n 

0.789(12) 

— 

2.307(94) 

— 

-0.228(773) 

— 


0.778(9) 

0.954(4) 

3.413(96) 

1.032(25) 

-0.630(915) 

-0.979(456) 


0.775(8) 

0.942(2) 

4.442(110) 

1.349(16) 

-0.280(852) 

-0.787(266) 

^ccc 

— 

0.937(2) 

— 

1.609(12) 

— 

-0.655(182) 



0123456789 10 11 

t. 


0.959 

0.952 

0.945 


43 „ 


H-1-1-1-1-h 


43 „ 


H-1-h 


H-1-1-h 


0.938 - Q 


_J_I_I_I_I_L 


0123456789 10 11 

t. 



FIG. 2. Strange (filled) and charm-quark (empty) contri¬ 
butions to the EO form factor at the lowest allowed three- 
momentum transfer (q^ =0.183 GeV^). The contributions are 
shown for single quark and normalised to unit charge. The fit 
regions are = [4, 7] for the charm sector and ti = [6,9] for 
the strange sector. 


tion between the electric charge radii of the spin-3/2 
baryons as, ((r|)n* + /2 = Com¬ 
paring ((r|)n* + /2 = -0.011(8), as obtained 

from such an estimation and the computed charge ra¬ 
dius of = —0.012(6), this relation seems to 

hold nicely. This implies that the individual quark con¬ 
tributions for each flavor to charge radii are similar for 
all baryons we consider here and their radii differ due to 


FIG. 3. Same as Fig. 2 but for the Ml form factor. Fit region 
is ti = [4, 7] for all cases. 


different quark compositions they have. 

It may be instructive to compare the contributions 
of the strange quark to the electric charge radii in the 
case of n/(ssc), ll*c(scc) and S*(ssu), E*(suu) baryons. 
The latter have been calculated in Ref. [5]. Such a com¬ 
parison would provide a better understanding as to how 
the charge radii are affected when the light quark is ex¬ 
changed by a charm quark. A comparison of s-quark elec¬ 
tric charge radii in 11/ - S* in Table V reveals the effect 
of changing the single u-quark by a c-quark: When the 
singly represented quark is heavier, the s-quark charge 
radius increases. In case of the 12/^ - E* baryons, the 
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FIG. 4. Same as Fig. 2 but for the E2 form factor. Fit region 
is ti = [4, 7] for all cases. Error bars are slightly shifted for 
clear view. 

TABLE IV. Electric charge radii of the fl, 12^1 ^*cc a-nd Uccc. 
Results are given in fm^. Quark sector contributions are for 
single quark and normalised to unit charge. Electric charge 
radii of spin-1/2 baryons are estimated through form factor 
fits as in Ref. [10]. Total electric charge radius of the spin-1/2 
He is estimated by the Eq. 24 since its electric form factor 
vanishes due to its zero electric charge. 



{r%)s 


(rE) 


[fml 

[fm""] 

[fml 

flc 

0.329(25) 

0.064(11) 

-0.177(18) 

^cc 

0.313(16) 

0.073(4) 

0.026(4) 

fl 

0.326(21) 

— 

-0.326(21) 


0.345(17) 

0.062(5) 

-0.189(12) 


0.348(16) 

0.078(3) 

-0.012(6) 

^ecc 

— 

0.084(3) 

0.168(5) 


doubly represented light quarks are changed to c-quarks. 
While the current precision does not allow a clear con¬ 
clusion, such a comparison again suggests an increase in 
the charge radius. 


Findings for the electric charge radii 

• Strange quark charge radii are insensitive to the 
baryon quark-flavor composition. 

• Charm-quark charge radii increase as the number 
of charm quarks increases in the composition of 
the baryon. 

• For singly charmed baryons, s- and c-quark charge 
radii are not affected by the spin-flip whereas the 
charge radii of doubly charmed baryons increase. 



FIG. 5. s-quark contribution to the electric charge radii of 
the spin-1/2 flc, flee and spin-3/2 fl, fl/ and baryons. 



FIG. 6. c-quark contribution to the electric charge radii of 
the spin-1/2 fie, flee and spin-3/2 fl/, fl/^ and fleec baryons. 


• We hnd the electric charge radius of Q to be 
{rl)n- = -0.326(21) fm^. 


D. Magnetic moments 

Magnetic moments of the baryons are related to the 
—0 value of the magnetic form factor Ml. We evalu¬ 
ate the magnetic dipole moment in units of nuclear mag¬ 
netons, 

Mb = Gmi(0)v -= Gmi(0)—(28) 

ZTTIq TTIq 

where mjv is the physical nucleon mass and is the 
baryon mass as obtained on the lattice. 

In order to make contact with the value of the mag¬ 
netic form factor at zero momentum transfer, G'mi(O), 
we apply the procedure in Ref. [26] and assume that the 




















































FIG. 7. Ratio of the quark contribution to the electric charge 
radii of the spin-1/2 fic, flee to spin-3/2 O, fl* baryons. qs /b* 
is a shorthand notation for the ratio, {r%Yg / where 

q is the quark flavor and B is the baryon. 



FIG. 8. The electric charge radii of the spin-1/2 flc, flee and 
spin-3/2 fl, Q* and fl/c baryons. Absolute values are shown 
for a better comparison. Data points denoted by a triangle 
indicate a negative value. 


momentum-transfer dependence of the multipole form 
factors is the same as the momentum dependence of the 
respective baryon’s charge form factor. For instance, the 
scaling of Gmi is given by 

G‘m‘.(I» = Gm,(9")^^. ( 29 ) 

where we consider the scaling of s and c quark sectors 
separately since each sector has a different scaling prop¬ 
erty. Gmi( 0) is then constructed via Eq.(24). 

Our numerical values for the magnetic moments are 
listed in Table VI, which are also illustrated in Figs. 9 and 
10. We find for spin-1/2 baryons that the contribution 
of a single quark to the magnetic moment significantly 
increases when it is doubly represented. The doubly- 
represented quarks make a spin-1 combination according 


TABLE V. Single strange quark contribution normalized to 
the electric charge radii of the 11/ and fl/c (normalized to 
unit charge) in comparison to that of the decuplet H* and 
S*. Decuplet values are taken from tables XL and XII. of 
Ref. [5]. 


Sn* 

Ss* 


5's* 

(rl) [fm"] 1 0.345(17) 

0.308(17) 

0.348(16) 

0.321(22) 


TABLE VI. Magnetic moments of the 17, 17/, and I7ccc. 
Results are given in units of nuclear magnetons, /ijv. Quark 
sector contributions are for single quark and normalised to 
unit charge. Magnetic moments of spin-1/2 baryons are esti¬ 
mated through form factor fits as in Ref. [10]. 




/^C 

d 


[mv] 

[MJv] 

[mjv] 

I7c 

0.979(47) 

-0.092(6) 

-0.688(31) 


-0.402(17) 

0.216(3) 

0.403(7) 

17 

1.533(55) 

— 

-1.533(55) 

17/ 

1.453(36) 

0.358(8) 

-0.730(23) 

H/c 

1.408(29) 

0.352(4) 

0.000(10) 

f^cce 

— 

0.338(2) 

0.676(5) 


to the Pauli principle and thus the dominant component 
is aligned with the total spin of the baryon, which leads 
to such an enhanced contribution compared to the anti¬ 
aligned singly-represented sector. A sign change is evi¬ 
dent due to the spin flip in case of the Oc and Qcc baryons. 
The s-quark contributions in the spin-3/2 sector have a 
slight tendency to decrease as the number of s-quarks 
decreases. In contrast to this, the c-quark contributions 
tend to decrease as the number of c-quarks increases. 

A comparison of the light sector to the heavy sector 
may reveal the differences between the systems having 
the same spin but different composition of quark flavors. 
For instance, we may compare the s-quark magnetic mo¬ 
ment in the charmed baryons to that in the light decuplet 
baryons, namely the 17/ (ssc) and S* (ssu) baryons or the 
I7/g(scc) and I]*(suu) baryons. Such a comparison would 
help us to understand the effects of changing a light quark 
by a charm quark. Magnetic moments of light decuplet 
baryons have been calculated in Ref. [5] with quenched 
lattice QCD. In Table VII, we quote the results obtained 
with the smallest pion mass, m,r = 300 MeV. As the 
work in Ref. [5] was done on quenched lattices at a much 
heavier pion mass, a quantitative comparison is not much 
useful. However, a qualitative comparison can be made 
as follows: The s-quark contributions to the magnetic 
moments of the charmed and light decuplet baryons are 
different. The charmed baryons have smaller magnetic 
moments than the light baryons. 

Similarly to the electric charge radii in Sec. Ill C, the 
effect of the quark spin configurations on the quark mag¬ 
netic moments can be further studied by considering the 
ratios of the quark-sector contributions to spin-1/2 and 
spin-3/2 baryons. The numerical values are given in Ta¬ 
ble VIII together with the octet-decuplet ratios extracted 
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FIG. 9. s-quark contributions to the magnetic moments of 
the spin-1/2 Qc, flee and spin-3/2 fl* and Qcc baryons. 
Absolute values are shown for a better comparison. Triangle 
symbol indicates a negative value. 



FIG. 10. c-quark contributions to the magnetic moments of 
the spin-1/2 flc, flee and spin-3/2 fl/, and fl^ee baryons. 
Absolute values are shown for a better comparison. Triangle 
symbol indicates a negative value. 


from Refs. [2, 5]. The ratios are also illustrated in Fig. 11. 
It is clearly seen that both the individual s-quark and c- 
quark contributions are enhanced in the case of spin-3/2 
baryons. The s-quark magnetic moment ratio in the case 
of doubly strange baryons, , and the c-quark mag¬ 

netic moment ratio in doubly charmed baryons, C^/ „„ , 
are consistent with each other. Such a behaviour is also 
observed in the case of the singly strange and the singly 
charmed baryon ratios, suggesting that the difference be¬ 
tween the spin-1/2 and the spin-3/2 baryons is almost 
independent of the quark flavour. 

We can further extend the comparison to octet and 
decuplet sector by including the E, S, E* and S* baryons’ 
strange sector magnetic moments extracted in Refs. [2, 5]. 
For the ease of discussion we compile the numerical values 


TABLE VII. Strange quark contributions to the magnetic mo¬ 
ments of the fl/, E* and E*. Decuplet baryon results 
are calculated in Ref. [5] on a quenched configuration with 
rrin = 300 MeV. All contributions are for a single strange 
quark of unit charge. 


Sn* 

Ss* 



fj, [hn] I 1.453(36) 

1.725(77) 

1.408(29) 

1.750(10) 


TABLE VIII. Ratios of the quark magnetic moment contri¬ 
butions in l/2^/3/2^. Octet/decuplet ratios are extracted 
from the numerical results available in the Refs. [2, 5]. All 
values are ratios of a single quark contribution of unit charge. 





Sh/h. 


0.674(34) 

0.615(10) 

0.703(50) 


'S'n„/n;^ 


'S'e/e* 


0.286(13) 

0.258(18) 

0.245(10) 


in Table VIII. It is seen that the s-quark contributions 
are similar for the doubly strange baryons and 

S'h/h*- However in the case of singly-strange baryons, for 
instance comparing the ratio to S's/s* ratio, we 

observe that the increase in the s-quark contribution (as 
the baryon spin changes from 1/2 to 3/2) is larger when 
the s-quark is accompanied by a uu component rather 
than a cc component. Interestingly, if we make a similar 
comparison for the singly represented quarks in 
and it is seen that the increase in single quarks 

contribution to the magnetic moment is less sensitive to 
the accompanying quarks flavor. 



FIG. 11. Ratios of the quark contributions to the magnetic 
moments of the spin-1/2 Qc, Dec to spin-3/2 D, D/ baryons. 
Rightmost blue data points are octet/decuplet ratios calcu¬ 
lated using the nriTr = 300 MeV quenched simulation results of 
Refs. [2, 5]. qB/B* is a shorthand notation for \i?bIP^b* where 
q is the quark flavor and B is the baryon. Absolute values 
are shown for a better comparison. Data points denoted by a 
triangle indicate a negative value. 

We calculate the total magnetic moments by combin- 
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ing the quark sectors via Eq. (24). The numerical values 
are given in the third column of Table VI and an illustra¬ 
tive comparison is made in Fig. 12. We find the magnetic 
moment of the baryon to be fiQ- = —1.533 ± 0.055 
fXN, which is smaller in magnitude than the experimental 
value, = —2.02 ± 0.05 /xat [25]. Magnetic moments 
are sensitive to the mass of the baryon. One of the rea¬ 
sons for this discrepancy can arise from the difference be¬ 
tween our mo = 1.790(17) and the experimental value, 
mo = 1.673(29), which is around 7%. Compared to 
the other lattice determinations that use the three-point 
function method, our value is slightly smaller (in magni¬ 
tude) than the quenched result, /ro- = —1-697 ± 0.065 
Hm, of Boinepalli et.al [5] and agrees with the Alexan- 
drou et.al’s extrapolated value, /ro- = —1-875 ± 0-399 
pLN [7] within one sigma error. In Ref. [27] magnetic mo¬ 
ment of O has determined to be /ro- = —1.93 ± 0.08, by 
a background field method on m^ = 366 MeV lattices. 

Magnetic moments of flc and il* are very close to each 
other suggesting that the spin flip of the charm quark has 
a small effect, as one would expect from a heavy-quark 
spin symmetry perspective. Based on a quark-model in¬ 
terpretation one would expect the magnetic moment of 
flc (flee) to be similar in magnitude to that of ft* 

While in the case of fie our finding is consistent with such 
an expectation, the magnetic moments of flee and fl*^ dif¬ 
fer drastically, the latter having a completely vanishing 
magnetic moment. The difference between the fie and 
fl* is that the c-quark is anti-aligned with the ss compo¬ 
nent in fie whereas it is aligned in fl*. Combined with 
their electric charges, quark sectors add constructively 
for fie and destructively for fl*. These two different be¬ 
haviours occur in such a balanced way that the magnetic 
moments of the fie and fl* end up to be similar. In case of 
the doubly-charmed flee and fl*^ however, the interplay 
between the electric charges and the number of quarks 
breaks the balance and lead to magnetic moments for 
flee and fl*j, that differ significantly. 


Findings for the magnetic moments 

• Quark sector contributions amongst the spin-3/2 
baryons are similar to each other, consonant with 
the quark-model expectations. 

• Magnetic moments of the strange and charm 
quarks in spin-3/2 charmed baryons are larger 
than spin-1/2 baryons having a similar quark- 
flavor composition. 

• Magnetic moment of the fl“ baryon is found to 
be, = —1-533 ± 0.055 /x^v- 

• flc, fl/ have similar magnetic moments in magni¬ 
tude. 



^GC ^ ^CC ^GCC 


FIG. 12. Total magnetic moments of the spin-1/2 flc, flee 
and spin-3/2 fl/, fl/^ and flcee baryons. Absolute values are 
shown for a better comparison. Data points denoted by a 
triangle indicate a negative value. 


• Magnetic moment of the fl/,. vanishes unlike the 

flcc- 

• As compared to the light decuplet sector, strange- 
quark contributions to the magnetic moments of 
spin-3/2 charmed baryons are smaller. 

E. Electric-quadrupole form factors 

The electric-quadrupole form factors of spin-3/2 
baryons provide information about the deviation of the 
baryon shape from spherical symmetry. In the Breit 
frame, the quadrupole form factor and the electric charge 
distribution are related as [28], 

QE2iy) = Ml J d^ry{r){iz^ - r^)x/>(r), (30) 

where 3z^ — r is the standard operator used for 
quadrupole moments. A positively charged baryon has 
a prolate (oblate) charge distribution when quadrupole 
form factor is positive (negative). 

As in the case of the EQ and Ml form factors, we 
estimate the E2 form factor by the plateau approach. 
We compute and extract the s- and c-quark sector con¬ 
tributions individually. E2 form factors in lattice units 
are shown in Fig. 4 and the numerical values, in units 
of e/M^, are given in Table IX. Unfortunately, low sig¬ 
nal/noise ratio does not allow us to conclude about fl 
and fl/ baryons. In the case of the heavier fl/^ and 
flccc baryons, however, the statistical precision is conclu¬ 
sive and it is possible to make a prediction about their 
shapes. fl/+ and fl++ have negative E2 moments thus 
their charge distributions deform to an oblate shape. 
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TABLE IX. E2{Q^) results for the Q, fl*, and Qccc at 
= 0.183 GeV^. Values are given in units of [e/m^]. Quark 
sector contributions are for single quark and normalised to 
unit charge. Last column is calculated by the Eq. 24. 



E2{Q^)s 

F2(g")c 

A2(g") 


[e/m^j 

[e/m1 

[e/m1 

n 

-0.337(1.142) 

— 

0.337(1.142) 


-0.371(539) 

-0.577(269) 

-0.137(352) 


-0.091(277) 

-0.255(87) 

-0.310(128) 

^ecc 

— 

-0.136(38) 

-0.273(76) 


Findings for the quadrupole form factors 
• and Q++ have oblate charge distribution. 


tent while keeping the m,r constant. Since this ap¬ 
proach is currently beyond our computational abil¬ 
ity, we resort to changing the while keeping 
the spatial extent constant. Rather than a thor¬ 
ough check, this approach is suitable to inspect the 
severity of the finite-size effects. We compare our 
Qc and Qcc results obtained on = 0.13781 lat¬ 
tices to the ones we have calculated previously on 
^u,d _ 0.13770 and 0.13754 lattices [10] for which 
ttIttL = 4.3 and niT^L = 6 respectively. Numerical 
results are given in Table X and shown in Fig. 13 for 
a clear comparison. Agreement within error bars 
suggests that the finite size effects are not severe 
for the charmed baryons although the = 2.3 
is well below the empirical value. 


F. Notes on systematics ^ 

V 

1. It is known that the Clover action has 0{amq) dis¬ 
cretisation errors and therefore the simulations may 
suffer from uncontrolled systematic errors when 
this action is employed for heavy quarks such as V 

the charm quark. In Ref. [9] we estimated the ef¬ 
fect of the discretisation errors arising from valence 
Clover charm quarks in a doubly-heavy Sec sys- 
tern to be small. Since we employ a similar ac¬ 
tion and formalism in the current work, we expect 
the discretisation errors to be small. However it 
is useful to check the significance of such system¬ 
atic errors once again with the current data set. 

We have repeated our simulations with parame¬ 
ter Kc = 0.1256 which leads to a decrease in the 
spin-3/2 charmed baryon masses by approximately /b 
3% MeV from those in Table II. Nevertheless such V 
a change in the charm-quark hopping parameter 
affects the EQ and Ml form factors by less than 
1% and ~ 3%, respectively. This change is prop¬ 
agated to the electric charge radii and magnetic 
moments as a change by less than 3%. As an illus¬ 
trative case, we quote the results for fl/ and Qcce 
baryons. We find (r|, Qccc)k,=o .1246 = 0.170(9) 
fm^ as compared to (r|, Occc)k,=o .1256 = 0.175(10) 
fm^ and = —0.696(50) as compared 

to = —0.712(50) ptN for 43 measure¬ 

ments. We note that this exercise gives only an 
indirect probe of discretization errors. A dedicated 
analysis with different lattice spacings is needed to 
quantify to which extent our results are prone to 
such errors. 



2.3 4.3 6 

ni_ L 
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1 1 
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FIG. 13. Electric charge radius and magnetic moment of fic 
and lice baryons on different m,rL values. Since the tic baryon 
has zero net electric charge we give the strange and charm 
quark charge radii. 


2. The set of gauge configurations we have used in this 
work has an = 2.3 value below the empiri¬ 
cal bound of 4, for which the finite volume effects 
are considered to be significant. In order to quan¬ 
tify the finite-volume effects, one should repeat the 
computation on lattices with different spatial ex- 


3. In order to extract the charge radii and magnetic 
moments of the spin-3/2 baryons, we use only two 
sets of momentum values, namely [Q^ = 0, = 1], 

and employ Eqs. (27) and (29). A more frequently 
used, alternative approach is to perform a full form- 
factor fit on higher Euclidean momentum region. 
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TABLE X. Electric charge radius and magnetic moment of 
fie and ficc baryons on different m,rL values. 


m^L 

2.3 

4.3 

6 

{r%)s [fml 

0.329(24) 

0.313(36) 

0.320(28) 

He {r%)c [fm^] 

0.064(11) 

0.061(10) 

0.076(13) 

M [fiiv] 

-0.688(31) 

-0.640(55) 

-0.621(44) 

Hcc {te) [fm"^] 

0.026(4) 

0.029(6) 

0.040(6) 

/i [^jv] 

0.403(7) 

0.402(15) 

0.400(11) 


The method we use allows us to estimate the ob¬ 
servables of spin-3/2 baryons more precisely since 
the higher momentum insertions introduce larger 
statistical errors. Although a fit to a set of 
values is desirable to estimate the observables, the 
values extracted by Eq. (27) or Eq. (29) agree with 
those obtained by a form factor fit. As an illus¬ 
trative example we can compare our values for the 
spin-1/2 Oc baryon given in Table XL For the spin- 
1/2 baryons we perform a dipole fit to their form 
factors as outlined in Ref. [10], which are quoted 
in the middle column. The values given in the 
rightmost column are extracted using Eq. (27) for 
comparison. Considering the current precision, we 
conclude that both approaches agree. 


TABLE XL Electric charge radius of Qc extracted by a dipole 
fit to the electric form factor and by use of Eq.(27). 



Dipole fit 

Eq.(27) 

[fml 

0.329(24) 

0.338(26) 

(4)c [fm^] 

0.064(11) 

0.067(11) 

(4) [fml 

-0.177(18) 

-0.180(20) 


4. Another source of systematic error is due to differ¬ 
ent fit strategies. One can either perform fits to 
the data sets of individual quark sectors and then 
combine the fit results via Eq. (24) to construct 
the baryon properties or, first, the data sets of the 
quark sectors can be combined to perform a fit to 
extract the baryon properties directly. In princi¬ 
ple both methods should lead to the same result, 
however a difference may arise due to statistical 
fluctuations. Such fluctuations cancel when data 
sets which are correlated with each other are com¬ 
bined, leading to a better estimated value especially 
if the value is close to zero like in the case of the 
When possible we choose to follow this 
strategy to perform our fits, e.g. spin-1/2 values 
quoted in Tables IV and VI. One caveat is that, in 
the case when the electric charge of the baryon is 
zero, combining data sets naturally leads to a zero 
electric charge radius. In such cases (e-g- (r|)nj 
we do not extract the observable from a combined 
data set. One can check whether two approaches 
agree with each other by comparing for example 
= —0.688(31) fj,N from Table. VI (fit to com¬ 


bination of quark sectors) with = —0.714(35) 
liN as obtained via Eq. (24) (fit to individual quark 
sectors). For the spin-3/2 sector, since we do not 
have more than two values, a form factor fit is 
not possible. Then we combine the fits to the data 
sets of individual quark sectors by using Eq. (24). 

5. When computing the matrix elements, one must 
optimize the temporal separation of the source and 
sink so as to minimize the excited-state contami¬ 
nation but also to obtain a good signal with mini¬ 
mum amount of statistics. For the electromagnetic 
form factors it has been shown that a separation 
of around 1 fm is sufficient [7, 10, 29, 30]. One of 
the strategies to extract less contaminated values 
is the so-called summed operator insertions (SOI) 
method [31]. In the usual plateau method a generic 
ratio has the form, 

R(t2, ti;p', p; r) = i?G + C> 

A'(t2-q)^ ^ 

where the Rq is the ground state value of the ratio 
and the first excited state contributions are sup¬ 
pressed proportional to ti and {t 2 — ti). A and 
A' are the energy gap between the ground and first 
excited state of the source and sink baryons respec¬ 
tively. In the SOI method, one sums the ratio in ti 
up to t 2 so that it assumes the form, 

*2 

R{t 2 , ti;p', p; r) = RgM + c(A, A') 
ii=o (32) 

+ 0(t2e-^‘^) +o(t2 e^'‘^) , 

where the c(A, A') is a constant. First excited-state 
contributions are now suppressed by t 2 ^ which is 
larger than ti or (t 2 — ti). It is possible to calculate 
the ratio with different source-sink separations and 
extract the ground state value Rq from the slope of 
a linear function in t 2 - Heavy-quark spin symmetry 
suggests that as the mass of the quarks increase, the 
energy gap between the ground and excited states 
of the baryons decrease. In Fig. 14 we show the EO 
form factors of and baryons, for which we 
expect the excited state contamination to be the 
severest. A comparison between the electric charge 
radius values extracted by plateau and SOI meth¬ 
ods are given in Table XII. Agreement between the 
results suggests that the excited-state contamina¬ 
tion is under control for t = 1.09 fm separation. 


IV. SUMMARY 

We have calculated the electromagnetic form factors 
of the H, H/, and flccc baryons at the lowest allowed 


+ o(e 
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" £2^^(strange sector) 


Form factor values 

extracted from SOI 




AeoWT 



• Q-2 - 0 

Pic 

Pccc 

0.769(59) 

0.936(22) 

0.934(15) 




- sector) 


secior/ 



• Q'2 - 0 

A Q-2 . 1 



• Q-2 - 0 

A Q-2 1 


10 12 14 10 12 14 

^2/^ ^2/^ 


FIG. 14. i?0 form factors of and IIccc baryons for two 
momentum insertions as obtained from summed operator in¬ 
sertions method. Results are for 43 measurements. Values 
given in the table are for = 1. 


to the number of the charm quarks that compose the 
baryon. Spin flip has a significant effect on the doubly- 
represented quark sectors so that the charge radii of the 
doubly-represented strange and charm quarks increase 
when the spin of the singly-represented quark is flipped. 

In case of the magnetic moments, quark-sector be¬ 
haviours change drastically between the spin-1/2 and 
spin-3/2 charmed baryons such that the quark-sector 
contributions to the magnetic moments of spin-3/2 
baryons get enhanced. Magnetic moments of flc 
n* baryons are found to be similar indicating a negligi¬ 
ble spin-flip effect by the singly-represented charm quark, 
n*,,, baryon has a vanishing magnetic moment unlike the 
spin-1/2 flee baryon. Strange-quark contributions to the 
magnetic moments decrease within the charmed baryons 
compared to the decuplet sector. 

Finally, we have been able to achieve a statistically 
significant data for the quadrupole moments of the 17/^ 
and Occe baryons to conclude that their electric charge 
distributions deform to oblate shape. 


TABLE XII. Electric charge radius of and flccc baryons 
extracted by the plateau and the SOI method. Results are 
compared for 43 measurements. 



Plateau 

* 

cc 

SOI 

Q 

Plateau 

zee 

SOI 

(r^)„ [fm^] 

0.332(43) 

0.359(113) 

— 

— 

(4)c [fm^j 

0.079(8) 

0.086(31) 

0.085(5) 

0.089(21) 

(4) [fm^j 

-0.005(16) 

-0.006(53) 

0.170(9) 

0.179(42) 


three-momentum value (q^ = 0.183 GeV^) on the lat¬ 
tices we use, and extracted their electric charge radii, 
magnetic moments and quadrupole moments. Based on 
the method outlined in our previous work [10], we have 
computed the electromagnetic form factors of the Oc and 
rice baryons also and extracted the electric charge radii 
and magnetic moments. For each observable we have 
identified the quark sector contributions as well as the 
baryon properties. 

We find that the electric charge radii of the strange 
sector is insensitive to the composition of the baryon, 
whereas the charm sector shows a slight dependence 
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